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Abstract 

Although an important issue in canonical quantization, the problem of representing the con- 
straint algebra in the loop representation of quantum gravity has received little attention. 
The only explicit computation was performed by Gambini, Garat and Pullin for a formal 
point-splitting regularization of the diffeomorphism and Hamiltonian constraints. It is shown 
that the calculation of the algebra simplifies considerably when the constraints are expressed 
not in terms of generic area derivatives but rather as the specific shift operators that reflect 
the geometric meaning of the constraints. 
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1 Introduction 



A key problem of canonical quantization is the incorporation of the classical constraints. 
In general, we have to expect that there are many different and in their outcome genuinely 
inequivalent ways to import classical symmetries into a quantum theory (for a recent example, 
see [H ) , and a famous example for the resolution of such ambiguities is the critical dimension 
of string theory. How to define and how to regulate the constraint operators is directly 
related to the representation of the constraint algebra in the quantum theory. Typically, the 
constraint operators are not uniquely determined by the classical constraints, and demanding 
well-definedness of the quantum constraint algebra is an important restriction on the choice 
of representation. 

Here we study the constraint algebra of vacuum general relativity in 3+1 dimensions 
in the framework of Dirac quantization, in which the classical constraints C are elevated 
to quantum operators C and imposed as operator relations on the wavefunctions, Cip = 0. 
Given a definition of the constraint operators, the question is whether there are anomalies 
in the constraint algebra, and if so, how to treat them. Currently there are no rigorous 
results about the constraint algebra of quantum gravity because of difficulties to define the 
constraint operators, e.g. issues of factor ordering and regularization arise. 

In the canonical formulation of vacuum general relativity there are two types of con- 
straints, the three-dimensional diffeomorphism constraint D and the Hamiltonian constraint 
H, which satisfy the following algebra: 

{D{v),D{w)} = D{C,w), (1) 
{D{v),H{N)} = H{C,N), (2) 
{H{M),H{N)} = D{g-'u;,), (3) 

where Ub = MdbN - NdbM, g"'' is the inverse three-metric, v and w are vector fields, and 
M and N are scalar densities of weight -1, all on a three manifold which we choose to be 
compact. 

The constraint algebra of general relativity has two important features. First, the algebra 
closes but only for structure functions containing one of the geometric variables in (^) . And 
second, neither D nor H form an ideal, and therefore one cannot construct a reduced phase 
space with respect to just one of the constraints. For example, the diffeomorphism constraint 
would form an ideal if the constraint on the right-side of @ was D, but as it stands, there 
are no invariant subspaces that are necessary for phase space reduction. 

For canonical quantization of general relativity a very fruitful approach has turned out 
to be the Rovelli-Smolin loop representation based on the Ashtekar variables [^, For 
a motivation of the loop representation and a discussion of its strong and weak points, see 
for example In this paper we focus on the technical problem of how to construct the 
quantum constraint operators in the loop representation, and present a formal calculation of 
the constraint algebra in the loop representation. 

The loop representation, in which states are functionals of loops, can be obtained from the 
connection representation, in which states are functionals of the Ashtekar connection, through 
the so-called loop transform. We introduce a for our purposes particularly convenient form of 
the constraints via the transform. The level of our discussion remains on the same heuristic 
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level as the original attempts to define the constraints ^, ^, ^. In particular, the loop 
transform is not defined rigorously. Although there now does exist a rigorous definition 
of the two representations and the loop transform based on distributional connections |^], 
which brings about certain changes to the whole formalism, it is not yet known how to 
treat the Hamiltonian constraint rigorously. Note that the constraint operators in the loop 
representation can also be obtained without the loop transform directly from the Wilson-loop 



variables, with the same result |10| 



The definition of the constraint operators requires a regularization, and we use a point- 
splitting regularization. The problem with point-splitting is that it introduces a background 
dependence which breaks diffeomorphism invariance and which survives in the limit that the 
regulators are removed. Further problems are that the constraint operators should act on a 
Hilbert space of states, but we do not have an inner product, and that the reality conditions of 
the Ashtekar formulation have not been implemented, and the constraint algebra of complex 
relativity might differ from that of real relativity. There have been important advances 



on these three issues (e.g. [11 1, [12|, and |]T3[, respectively), but for the above reasons all 
calculations of the constraint algebra must still be called formal. What can be explored, 
however, is what the current framework predicts for the constraint algebra, and in the loop 
representation this is not a trivial task. 

One can distinguish four approaches to the constraint algebra in the loop representation: 

1. Recall that in the connection representation, the constraint algebra, formally regulated 
by point-splitting, does close and one can hope that this feature is preserved by the 
loop transform (which is as we already pointed out still problematic for the Hamiltonian 
constraint). This justifies the route most commonly taken until quite recently, namely 
to postpone the treatment of the constraint algebra in the loop representation. Given 
the definition of the constraint operators in the loop representation, one proceeds to 
study their kernel, to look for an inner product, defines observables, etc. without taking 
the constraint algebra into account. 

2. Given the definition of the constraint operators in the loop representation, the constraint 



algebra is computed explicitly |14]. The technical difficulties encountered explain strat- 
egy 1. One problem is that the differential operator on the space of loop functionals 
that appears in the Hamiltonian constraint [^, ^, the area derivative, is not very well 



studied, but now the machinery is available ||T5[ that makes |14| possible. 



One attempts a two-step procedure, first solving the diffeomorphism constraint, and 
then defining the Hamiltonian constraint operator only on the diffeomorphism invariant 
states. This is how sometimes the original Rovelli-Smolin loop representation [Q] is 
interpreted, namely as a theory of operators on knot invariants. For example, in the 
rigorous approach based on measures on the space of connections modulo gauge [ p!^ , 
the measures are constructed to be diffeomorphism invariant first, and one later 
attempts to construct a "diffeomorphism invariant" Hamiltonian constraint operator. 
An analogous approach is taken in Loll's lattice gravity, which has led to important 
insights for certain geometric operators |17]. 

From the perspective of the quantum constraint algebra corresponding to (|5)-(^), we 
have to observe that one cannot simply first impose the diffeomorphism constraint 
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without consequences for the Hamiltonian constraint, since as mentioned above the 
constraints are intertwined such that there are no ideals. Let us give a simple-minded 
example. Suppose that D{v)iIj = and that H{N) does not map states out of the kernel 
of D{v), i.e. D{v)H{N)i; = 0. Then [D{v),H{N)] = 0, and by (|), H{C^N) = 0, 
which since it holds for all v and N implies H{N) = That is, if we want to 

work exclusively on the space of solutions to the diffeomorphism constraint, we have to 
consider the subspace that is invariant under the Hamiltonian constraint, and then the 
constraint algebra implies that this subspace must be the simultaneous kernel of the 
constraints. 

If we want to define the symmetries of quantum gravity with the help of the opera- 
tor constraint algebra, we therefore have to implement the constraint algebra before 
imposing the diffeomorphism constraint. 

In this context it is perhaps worth pointing out that one of the main attractive features 
of the loop representation, namely that a few non-trivial (formal) solutions to both 
constraints are known, is not derived from a two-step implementation of the constraints. 
Rather, one attempts to find the intersection of the kernels of the constraints on the 
unreduced space of loop functionals |19]. 



4. One introduces (almost trivial) matter variables that allow the definition of a preferred 
time-slicing. The Hamiltonian constraint operator then turns into a proper Hamilto- 
nian operator, and the Wheeler-DeWitt equation into a proper Schrodinger equation 
[ prfl . The main advantage is that the regularized Hamiltonian operator is diffeomor- 
phism invariant, which opens up a whole range of interesting topics in diffeomorphism 
invariant dynamics. The problem of representing the constraint algebra is easily solved 
because only the diffeomorphism algebra remains to be represented, the information 
about the Hamiltonian constraint is now contained in the Schrodinger equation. In this 
sense, a preferred time-slicing also defines a two-step procedure. One should remem- 
ber, and emphasize, however, that the construction of matter clocks works only locally. 
The resulting theory is therefore only an approximation to what is usually called full 
quantum gravity. It would be very interesting to gain some control over the approxima- 
tion, for example in reduced models. Even so, the prospect of diffeomorphism invariant 
dynamics is very interesting, and may be physically quite relevant. 

To summarize, one can either try to implement the classical constraints directly, 1 and 
2, or to give a special treatment to the Hamiltonian constraint, 3 and 4. For the latter, 
there are good motivations, like the observation that while the diffeomorphism constraint is 
linear in the momenta and therefore generates a type of gauge symmetry, the Hamiltonian 
constraint is quadratic in the momenta and does not possess simple gauge orbits. Also, to 
obtain the conventional interpretation of time in quantum gravity a special treatment of the 
Hamiltonian constraint may be necessary. 

Here we do not choose to treat the constraints differently, but want to explore as in 2 
whether the classical constraint algebra has a complete representation in the loop represen- 
tation of canonical quantum gravity. To state the outcome of the calculation, the resulting 
algebra maintains the structure of the classical algebra with a particular choice of factor 
ordering for D{g°'^u}})), which becomes necessary because of the presence of the metric. The 
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level of rigor is equivalent to that in the connection representation, and upon removal of the 
point-splitting there are no anomalies. Let us point out that the formal nature of the point- 
splitting regularization does not allow us to decide whether there actually are anomalies in 
the constraint algebra of quantum gravity. 

Our contribution is to show that starting with a different form of the constraints the 
result of Ijl^ can be obtained in a simpler manner. The simplification comes about by the 
observation that the constraints are not just based on generic area derivatives, but rather 
on more specialized geometric operators, certain shift operators. This makes the algebra 
managable to the extent that the point-splitting regularization can now be studied further 
along the lines of where the removal of the regulators is examined in detail for the 
constraint algebra in the connection representation. In terms of shift operators, the constraint 
operators may in fact be compatible with the rigorous regularization techniques coming from 
distributional connections, with which it is particularly hard to represent the field strength 
of the Ashtekar connection that directly corresponds to the area derivative via the loop 
transform. 

The paper is organized as follows. In section 2, we define the constraints in the connection 
and the loop representation and discuss the point-splitting regularization. In section 3, we 
introduce the basic loop derivative commutators. Section 4 contains the calculation of the 
constraint algebra in the loop representation. In section 5, we conclude with a few comments. 



2 The constraints and their representation 

The constraints of quantum gravity in the loop representation have been derived in various 
ways [§, |6|, 0, with essentially the same result |1C]. Since our method to compute the 
constraint algebra depends crucially on the form of the constraint operators, let us give a 
brief derivation via the loop transform. 

The Ashtekar variables are a connection j4^(x) and a vector density of weight one, 

both complex, on a compact three- manifold S. Tangent space indices are denoted by a, 6, . . ., 
internal indices are denoted by The internal gauge group is SU{2), and following |^] 

we choose generators r* such that 

[r\r-''] =e^^'V. (4) 

The algebra-valued variables are obtained by contraction, e.g. Aa = ^gj^ ■ The inverse metric 
is given by 

where g is the determinant of gab (insuring the correct density weight). 

The constraints of general relativity which satisfy the algebra (|l|)-(^) are 



B(v) = jd'xv'^E''Fl, + G(v\ (6) 
H(N) = IS'TNi^^E'^^E^^Fl^. (7) 



The vector constraint / S'x v"'E'^^F^f^ generates diffeomorphisms up to an internal gauge trans- 
formation, which is compensated by the term G{v), and a gauge depending term also appears 
in (^. (The sign convention for the Poisson brackets in (|l])-(|^ is opposite to that of ([§]).) 
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In the connection representation, wave functions are functionals of the Ashtekar connec- 
tion, ip[A], and the operators corresponding to the connection and the triad are represented 

by 

Ai{x) = Ai{x), ^»(x) = -^^, (8) 

where h = 1 and a complex i has been absorbed in the definition of i?"*. 

The first non-trivial issue we have to face is regularization. While the type of point- 
splitting that we use has been commonly applied in many places, let us proceed slowly since 
there are different prescriptions for in which order the various regulators have to be removed. 

The necessity for regularization arises at this point because the metric and the constraints 
are products of operators at the same point. We introduce a point-splitting based on a 
background metric and a regulator f^{x,y) satisfying 

limUx,y)=6Hx,y). (9) 

For the calculations that follow we fix a coordinate system and require that fe{x,y) is a 
smooth function of x — y. 

As discussed in |21, which refer explicitly to the constraint algebra, a 'full' point- 
splitting should be applied, that is, all points that appear in an operator product should be 
split. Also see [^^, where it is shown that only in a particular factor-ordering the connec- 
tion representation is related to the loop representation, and that the vector constraint in 
this factor ordering only gives rise to diffeomorphisms in the connection representation if a 
symmetric point-splitting, f^{x,y) = f^{y,x) is employed. Following these considerations we 
define the operators 

,fAx) = (10) 

DM = -Jd'x Jd%Mx,y)v-{x)^^^F:,{y) + G{v), (11) 
H.AN) = Jd^x jdS |d3z/,(x,y)A,(x,z)A^(x)6*^'=^^^^F,Vx). (12) 

For the particular calculations that follow, some of the regulators can be removed since 
there do not arise singularities that require them. (As just recalled, this is not true for the 
vector constraint algebra in the connection representation.) 

The transition to the loop representation is made via a formal transform, the loop trans- 
form 

V'M = jvAtiU^^[Al (13) 
where C/^ is the holonomy matrix of around the loop 77, 

Ur, = Pexp j ds f]''{s)Aa{ri{s)). (14) 

In order to transfer the operators that we are interested in in the connection repre- 
sentation, Oci to operators in the loop representation, O^, one performs a formal partial 
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integration for any occurrence of a functional derivative with respect to so that 

Ol^W = JvAtTUr^Ocip[A]= JvA{0^,tTUr,)^p[A]. (15) 

All that is needed for an explicit definition of as loop operator is a transfer relation of the 
operators on the Wilson loops that expresses the operation on the connection dependence of 
the Wilson loop purely as an operation on its loop dependence, 

O^tra^ = OitiUr,. (16) 

This construction is directly applicable to the metric and the constraints. The Wilson 
loop satisfies 

jJj^tvUr, = Jds5Hx,v{s)W{sMUssn), (17) 

^^trC/, = ^\s)FMs)MU,sTi), (18) 

where Ugs denotes the holonomy from r]{s) once around the loop. Variations of the Wilson 
loop with respect to the connection and the loop are not completely unrelated, and for the 
operators under considerations there exist transfer relations precisely because of that. In fact, 
(15) is the one transfer relation we need, relating the loop derivative on the loop representation 
side with the multiplication and insertion of a field strength on the connection representation 
side. 

At this point it turns out to be quite advantageous to introduce a new piece of notation. 
Let Tg be the loop operator that inserts a generator at parameter value s into the holonomy, 

TitiU = tiUssT\ (19) 

Partial integration as in (|^) and evaluating the functional derivative with respect to the 
connection as in leads to 

DivMr]] = lvA(^ldsv-ivism\s)FMs)mtrU,^m] (20) 
H,e{N)i^[v] = JvA (^Jd^x Jds Jdtf,{x,ij{s))Mx,i]{t))N{x)r{s)ii\t) 

e^^'T:T^F^,{x)trU,) m (21) 
gee'i^mv] = fvA( fds [ dtUx,ii{s))fAx,v{m''{s)v\mn^^U^ij[A]. {22) 



We have removed the regulator in D{v), and G{v) does not contribute since the Wilson loops 
are gauge invariant. 

A single insertion operator could be transferred to the loop representation by introducing 
functionals of loops with a marked point, Tg^[ri]. Such an extension is not necessary if the 
insertion operator is combined with a field strength as in D(v), or if there are two insertions 
as in g"''', or if there are two insertions and a field strength as in H{N). The double insertions 
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that arise combine due to the trace identities for SL{2, C) matrices to the following rerouting 
operations, 

rXiiU = ^tiU -UiUtsirUst (23) 
e'^^TlTiiiU = ]^{iTUstirUtsr^ -iTUstThiUts) (24) 

where Ugt denotes the parallel transport from s to t going around the loop in the positive 
direction. That is, if s > t, Ust = UsiUot- Note that the resulting holonomies in ( ^J2^ ), e.g. 
trUst, in general refer to open paths, which however reduce to closed loops in the limit that 
the regulator is removed. 

As demonstrated, taking two derivatives with respect to A^^ leads to the well-known 
rerouting of the loop at intersections in the loop representation. The rerouting in g"'^ and 
H{N) is by no means a trivial side-effect but actually crucial for the definition of the operators 
and the closure of the constraint algebra. The Hamiltonian constraint includes differentiation 
as well as rerouting. The trace identities capture the fact that the internal gauge group is 
complexified SU (2) and not some other group. 

The rerouting has to be denoted in some way, and we find it simplest to not resolve the 
insertion operators into rerouted loops whenever possible. For example it is then irrelevant 
whether s < t ov t < s, and where the parameter origin of the loop is (as opposed to |l^]). 
This turns out to be an important technical point, since the use of insertion operators allows 
us to separate the rerouting operations from the other calculations and lead to a significant 
simplification. 

Whether we work with reroutings or insertion operators, we have to deal with the case 
s = t. Note that for a given s, tvUtsT^UstT^ is a function in t with a finite step discontinuity 
at t = s. The integrations in the metric and the Hamiltonian constraint are evaluated by 
considering left- and right-sided limits, Jds Jdt = Jds Jq dt + Jds /j+ dt. 

Note also that the product T^T^tvU has no intrinsic meaning, since it is not clear where the 
second insertion has to take place, before or after the first, but the left- and right-sided limits 
TgT^± are well-defined. Therefore, we define insertion operators on general loop functionals 
that produce functionals of loops with a marked point, Tgijj[ri], and that inherit certain trace 
identities from the r*. In particular, 

[T;,r/] = fors/t, (25) 
[Ti,rf] := Tl.Ti^-Ti+Ti_ (26) 
[TlTi] = e'^^Tl (27) 

The transfer relation ([T^ can be directly applied for D{v) and g'^^,{x). For H{N) there 
remains the problem that, as given in (^), the field strength is evaluated at x and not on 
the loop as required in (|l8|). One can either introduce at this point a more general type of 
loop derivative, namely a path dependent area derivative (cmp. section 3.1), or use that in 
the integrand 

F:,(x) = FMt)) + 0{e') ^ FMt)). (28) 

All our calculations (and those of p!^) are performed only to leading order in the point- 
splitting. 
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The transfer relation (|T^, the commutator (|2^), and attaching the field strength to the 
loop (|2^), allow us to arrive at our final form for the metric and the constraints, 

D{v) = Jdsv^rjis))^^, (29) 

Hee'iN) = jSx jds JdtUx,7j{s))fAx,v{t))r{s)N{7^{tm[^^^^^^^^^ (30) 

g${x) = Jds jdtUx,ii{s))fAx,v{tW{s)fl\t)TiTi, (31) 

where in (^0|) we have assigned the marked point property also to the loop derivative. 

The constraints in the loop representation that we have derived are equivalent to the 
standard result 0, Apart from notational differences for the rerouting, and order of 
e' differences in the regulators, there is, however, one very important technical difference. 
The loop derivative that we introduce for the transfer is not the area derivative, but just a 
special case thereof, the ordinary functional derivative with respect to the loop, which has the 
geometric interpretation of an infinitesimal shift operator. Considering the other approaches, 
it is not clear why the functional derivative should suffice, and we show below that for the 
removal of the regulator one is forced to introduce the area derivative at generic kinks of the 
loop. In order to arrive at the above form we had to take special care that the tangent vector 
to the loop is always at the same point as the field strength as in (|l^). 



3 Loop derivatives and some basic commutation relations 

Before moving on to the computation of the commutators of the constraints, we take a more 
detailed look at the loop derivatives that appear in the constraints, and analyze the basic 
commutators of the derivatives and the rerouting operators. 



3.1 Commutators of loop derivatives for smooth loops 

The basic commutator for the computation of the constraint algebra is that of two functional 
loop derivatives, 

^W{s)'Wit)^^^' ^^^^ 
While this commutator vanishes, the commutator of two area derivatives does not vanish, 
which is a good reason to attempt to rewrite the Hamiltonian constraint in terms of functional 
derivatives. Let us discuss this important point in more detail. 

The area derivative of a loop functional [rf\ is defined by appending to r/ an infinitesimal 
loop 7'' with area element cj"^(7'') = 0(l/(5^). In its general form, the area derivative depends 
on a path (and its inverse 7r°) from the point o at which all loops r/ are supposed to be based 
to the point x where 7*^ is attached. The definition of the path dependent area derivative is 

Aaf,(vrJ^[r?] = hm , (33) 
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Figure 1: Parameter dependent area derivatives that act on different points of the main loop 
commute since it does not matter in which order the infinitesimal loops are inserted. 

where juxtaposition of loops denotes attachment of the loops at the base point. From that 
definition follows the basic commutator used in [ p!^ , 

[^ab{<), ^cdi^l)] = ^ab{K)[^cd{^l)l (34) 

where the brackets on the right-hand side indicate action on the path dependence in vr^ 
only and not on the loop functionals. Note the peculiar form where the commutator of two 
derivatives is the derivative of a derivative. 

The area derivative that typically appears in the derivation of the Hamiltonian constraint 
does not depend on arbitrary paths but on portions of the loop on which the area derivative 
acts. From definition ( |3^ ) with tTq = rj^, x = r]{s) we have a definition for the parameter 
dependent area derivative, 

A.,(,)«„) := AMrPyym = lim /'^';:,^!-'"'''! . (35) 

As remarked in the parameter dependent area derivative is naturally more restricted 
in its applicability than the generic, path dependent area derivative, but let us point out 
that nevertheless it is all we need for the commutators of the constraints. To be more 
specific, explicit calculation of the infinitesimal loop variations show that [Aab{s), Acdit)] is 
not expressible in terms of a parameter dependent area derivative, but we also find that 

v\s)v''it)[Aabis),A,a{t)] = -d{s,t)-^^Aac{s). (36) 

One can peel of one of the tangent vectors obtaining a so-called covariant loop derivative of 
an area derivative, but removing both tangent vectors does not leave a single area derivative. 

In pictures, parameter dependent area derivatives commute if the infinitesimal loops are 
inserted at different parameters of the main loop, fig. 1, but there is a nontrivial contribution 
if one of the small loops is inserted onto the other, fig. 2. 

To motivate the relation between the parameter dependent area derivative and the ordi- 
nary functional derivative with respect to a loop, compare the transfer relation (^) with the 
Mandelstam equation, 

Aab{s)tTUr, = rMs))TltlUr,. (37) 

All that is missing is the contraction with a tangent vector, and one can show that for smooth 
loops, 

-^ = ^^(.)A„,(.). (38) 
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Figure 2: Two area derivatives at the same point do not commute since there is a non- 
vanishing contribution when one of the small loops is inserted onto the other. 



And as we already pointed out, the commutator of two functional derivatives vanishes. 

To summarize, the three types of derivatives that we consider are in increasing order of 
generality the functional loop derivative, the parameter dependent area derivative, and the 
path dependent area derivative. They are related by (^) and (|38|). The basic commutators 
are the more complicated the more general the derivative is, compare (^), (^), and (|3^. 

That suggests that using parameter dependent area derivatives already offer a simplifi- 
cation over path dependent area derivatives, but using functional loop derivatives is even 
simpler. In section 4, we compute the commutators of the constraints. It turns out that most 
of the calculations can be performed on the level of the functional loop derivatives, but at one 
point we fall back onto parameter dependent area derivatives. While it is not proven that it 
is not possible to perform the calculation exclusively in terms of functional derivatives, it is 
a convenient approach. 

As already emphasized, however, the Hamiltonian constraint is not just a derivative op- 
erator, but contains a rerouting as an additional complication, which we address in section 
3.2. 



3.2 Commutators of loop derivatives on piecewise smooth loops 

The loop representation allows continuous, piecewise smooth loops, which means there may 
be kinks, f]'^{s~) / ?7"(s"*"). The area derivative is obviously well-defined at kinks, e.g. ( p7| ) 
(that is, kinks do not prevent a loop functional from being area-differentiable). As evident 
from jf^ij^ = r]^{s)Aab{s), the functional loop derivative is ill-defined at kinks where "q^is) 
does not exist. 

For loops with kinks, integration around the loop as in Jdsfi"'{s) is defined in terms of 
left- and right-handed limits the same way we resolved the ambiguity in T^T^ for s = t. 
Even if one assumes that the loop argument of iplr]] is smooth, action by the metric or the 
Hamiltonian constraint introduces kinks at intersections. Recall that 

Tn ■= T-L-.Ti - T-^Ti, (39) 



and hence the loop derivative is kept away from the kinks introduced by the rerouting. 



11 



But notice that there are now two hmits involved. Suppose that there is a kink at Iq. In 
the limit that the regulators are removed, one of the terms in the Hamiltonian arises for t in 
( |30| ) close to to- We impose that the limit in the has to be taken before the limit of t ^ to- 
That means that always either t~ < t < < to or to < t~ < t < t"*". This definition of 
the integrals in the metric and the Hamiltonian constraint together with the point-splitting 
regularization, which led to our prescription for the positioning of the loop derivative, insure 
that the constraints are unambiguously defined for kinks. 

We therefore have for the commutators between reroutings and functional loop derivatives 
for integrations of the type that appears in the constraints (i.e. with the appropriate left- and 
right-sided splits in the range of integration), 

ds Jdtf{s,t)[Ti,Tn=0, (40) 

/../dt/.n,(.,t,n)[^,7^[^,I?]] = 0, (41) 

ds Jdt jdu Jdvh{s,t,u,v)[T^[j^yTn,n[^-^ (42) 

for continuous functions /, g, and h that maintain reparametrization invariance. For example, 
( |40| ) follows from [r*,T/] = Sgt^^-'^Tg as long as f(s,t) is assumed to be continuous. That is, 
there are contributions from the commutators defined by the left- and right-sided limits, 
but the integral is zero as long as these contributions are finite and have support only on 
sets of measure zero. Of course, the reason why the commutators of the constraints are not 
trivially zero because of (^0|)~(|4^) is that functional differentiation can lead to distributional 
coefficients /, g or h. 

Now we are ready to compute the constraint algebra. In section 2 we promised a demon- 
stration that the unusual form of the Hamiltonian constraint, that it involves only functional 
derivatives, reduces to the standard form with area derivatives at kinks when the regulator 
is removed. We do not actually remove the regulator in this paper, since this would involve 
several case distinctions for smooth portions of the loop, kinks and intersections, but as an 
important example consider the Hamiltonian constraint for a loop with intersections. Sup- 
pose r/(so) = ??(io) for So < to. One of the terms that arises due to the rerouting refers to 
V'[f?soto], i-e. for a generic intersection there is a kink at so, to. First notice that Aafe(so)^[?/soto] 
is well defined, in fact, 

Aa6(so)V'[??soto] = ^ab{tQ)^[risoto\- (43) 

Our definition of the Hamiltonian constraint gives rise to four terms with a functional deriva- 
tive given by the four orderings of to, t, and t^. For example, for to < < t < t+, 

t>t^to Jij^^it^^^'^'otoritot-Vt-tl = r]\4)^ab{to)ip[r]soto]- (44) 

The point is that in the limit that the regulator is removed we recover the standard result (e.g. 
Q) that in the generic case, when the loops do not run smoothly through the intersection 
but the tangentvectors ?7"(so ), rj'^isQ), ?7°(to ), and 77" (tg") are all different, there are terms 
for which f]"" {s)r]^ (t) Aab{t) cannot be combined into loop derivatives since the area derivative 
and tangent vectors are located on different legs of the intersection. 
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The form of the Hamiltonian constraint ( pO| ) , and the prescription for the order in which 
to take the two limits related to reroutings is adopted precisely because it allows us to always 
keep the area derivative and one tangent vector on the same leg so that they combine to a 
functional loop derivative. 

4 The constraint algebra in the loop representation 

4.1 Commutator of D(v) with D(w) 
Given (|32|), we immediately have that 
[Div),D{wmr^] 

(45) 

ds {v-{r^{s))daw'{v{s)) - w-{rj{s))dav'{v{s))) ^^^m (46) 

= DiC^w). (47) 

The same calculation can be performed after replacing 5/5ri{s) by if {s) /S.ah{s) , where the non- 
vanishing commutator of the area derivatives cancels the contribution from the variation of 
the tangent vectors. Much more involved is the direct use of path dependent area derivatives 



1 14], where the non- vanishing basic commutator cancels together with the variation of the 
tangent vectors and the path dependence in the area derivative when the Bianchi identity 
for area derivatives is used. The basic structure in all three cases is that the variation of the 
smearing vector fields gives rise to the relevant term in the commutator without contribution 
from the other variations. 

4.2 Commutator of D{v) with H,,i{N) 
By definition of the constraints and ( pl| ) we have 

[D{v),H,,,{N)] 

d'x jds jdt jduN{x)v^viu))-4^ {Ux,7jis))fAx,r]itm\s: 



57]"- (u) 

5_ 

5rf{t) ' 

d^x Ids jdt I duN{x)Ux,rj{s))fAx,7]it))ris) 



The functional differentiation of the tangent vector gives 

= ^^'(^(^)) = ns)dav''{v{s)), (49) 
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while for the vector field, 

Since in the limit that the regulators are removed r/(s) ~ these two terms cancel in the 
commutator, 

[D{v),H,AN)] 

Sx jds jdt JduN{x)v''{r,{u))^^^{Ux,v{s))fAx,vmv\s) 

For the functional differentiation of the regulators we obtain 

f S 
d'x j duN{x)v%7^{u))^^^ {Mx,^{s))fAx,vm 

d^x N{x){-v-{r]{s))fAx, Vit))dafe{x, vis)) - v''Ht))Ux, v{s))dafAx, vm (52) 

d^x {v''{v{s))Mx, r]{s))da{N{x)fAx, vm + v''{v{t))fe'{x, v{t))da{N{x)Mx, v{s)))) 

(53) 

d^x {v^{x)Ux, v{s))da{N{x)fAx, vm - da{v''{x)fAx, v{t)))N{x)Mx, vis))) (54) 

d'^x{v''{x)daN{x)-N{x)dav''ix))f,ix,rj{s))fAx,v{t)). (55) 

We performed two partial integrations in order to remove the derivative from the regulators, 
and we moved the vector field from the loop to x. Note that in general v'^{y)da6^{x — y) 
v'^{x)da6^{x — y), and that we first removed the partial derivative from the regulator before 
concluding v"'{ri(s)) ~ f"(x). 

Since is a scalar density of weight —1, CyN = v°'daN — Ndav"", and inserting (^) into 
dsip we obtain the operator version of the classical Poisson bracket, 

[D{v),H,AN)]^H,AJ^vN). (56) 

4.3 Commutator of Hs5'{M) with H,,,{N) 
By definition of the constraints and (^) we have 

[Hss'{M),H,AN)] 
>x fd% fds fdt fdu fdvM{x)N{y)fs{x,v{s))f5'{x,v{t))v''{s)Ti 



{UyMu))fAyMv))v'{u)) ,T/]r„^[^,r,^l 

-{M ^ N,5 ^ e,6' ^ e'). (57) 
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Functional differentiation of the tangent vector gives 

duMy,v{u))Ttj^^f,\u) = -±(f^{y,r^{t)m')6', (58) 

= mdcUy,vitm'5l-UyMt))^T,'6l (59) 
The second term does not contribute to the commutator since 

M{x)N{y)fs{x, r]{s))fs> {x, v{t))My, rjit))fAy, viv)) 

-N{x)M{y)Mx,rj{s))fAx,r){t))fs{y,v{t))f5'{y,v{v)) ^ 0. (60) 
Therefore, the functional differentiation of the tangent vector in the commutator gives 



d-'x jd-'y Ids jdt I duM{x)N{y)fs{x,r,{s))fs:ix,r,{t))d^Uy,r,{t))fAy,v{u)) 

6_ 

5rf>{u) 

{M ^ N,6 ^ e,6' ^ e ). (61) 



r?^(.)r)"(t)ri[r^r/][-^,r„^] 



To remove the partial derivative from the regulator, we perform a partial integration in y, 
and with approximations similar to (SO), 

>x J d'yM{x)N{y)fs{x,7]{s))fs'{x,i^{t))dafe{y,7l{t))fe'{y,v{u)) 

J dSN{x)M{y)Mx,r,{s))fAx,v{t))daf5{y,v{t))f5'{y,v{u)) 
fd^yu;a{x)fs{x,v{s))f5'{x,v{t))fe{y,v{t))fe'{y,v{u)) (62) 



where 



~ -uJa{viu))fs{ri{u),r]{s))f,{r]{u),r]{t)), (63) 



LOa{x) = M{x)daN{x) - N{x)daM{x). (64) 



{uja is a covector density of weight -2, and uJafsfe has weight 0.) In (^) we integrate over 
x and y keeping terms to leading order. Note that already in the definition (3C) of the 
Hamiltonian constraint one of the regulators can be removed without any renormalization, 
say f^i{x,r]{t)) S^{x,r]{t)). But after integrating over x, it is then not quite clear how to 
perform the necessary partial integrations in the commutator algebra. 

Functional differentiation of the regulators in the commutator (|57|) leads directly to par- 
tial derivatives of regulators, which can be removed analogously, and we obtain for the full 
commutator (after appropriate renaming of integration variables and contracted indices) 

[Hss'{M),H,AN)] 

'ds fdt fduu;air]{u))fsir^{u),^{s))fMu),v{tW{s)v''{t) 
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A priori it is not obvious how the right-hand side of the classical Poisson bracket D{g"'^LOi,), 
should be represented in the loop representation because of the operator product of the metric 
with the generator of diffeomorphisms. Comparing the above with the loop operators for the 
diffeomorphism constraint ( p9| ) and the metric (pT]), a natural guess is that the reroutings in 
( |65| ) combine to T/T/ g^b(^^-^ ■ This is indeed the case. 
The first two terms in (l65|) cancel since with (12^ 



' ^ e'^^TlTiT^ + e'^^TiTlT^] (66) 



= 0. (67) 

Note that the cancellation does not depend on whether j^^j^ is located at u~ or . 

However, the one remaining term cannot be simplified using the SU{2) commutator for the 
insertion operators, because resolving the order of the insertions puts the derivative between 
the insertion operators, 

AAA A A 

rr " rpkl rpjl _ " rpkrpj rpk " rpj _ T^J " rpk , rpj rpk " 

^^d7?^'(u)'^J' ''^ ''5i]\u+y^++ 5r]\u-) '' '' '' 5i]\u+y 

(68) 

Of course, there are other identities to work with, but let us first consider the situation in 
the connection representation, i.e. ■0[?7] = trC/^ and ^^/^^^ = ?7''(n)F^^(r/(n))T*. Then the 
rerouting simplifies according to 

The important observation is that this result cannot be transformed back to a functional 
derivative and some reroutings in the loop representation because in the first term Fl^^{r]{u)) 
is contracted into an insertion at s, not u. In the limit that the regulators are removed in 
(|65|), we have ri{u) ~ ??(s), and 

= ^^(.)A,.(.)T/r^ - TMj^y (70) 

As already discussed in section 3.2, the reroutings in the Hamiltonian may split a functional 
loop derivative into a tangent vector and an area derivative at different parameters. 

While by definition of the loop representation its states have to satisfy the loop transforms 
of all the SU{2) identities of the connection representation, i.e. in particular (|70|), let us also 
indicate how (^) can be derived using the SU{2) identities in the loop representation without 
resorting to the loop transform. For definiteness, consider the case s < t < u. The rerouting 
operations that appear in ( [70| ) are 

TiT/V'M = \m-li^[r]st^r]ts], (71) 
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Marking the area derivative by the insertion of a smah loop 7, the claimed relation ( [70| ) 
becomes 



tpiVst U -qus U r/t„7] - ip[r]tu U 7]us U 7]sa] + 'ip[r]us'ntu'nsa] - ipivtuVstilusl] 
= 2^[r]st U r]usi]tu7] - '^ipivtu U ??st??ns7] + V'[%t'7t«^ns7] - tpi'Husilst'ntul]- (73) 

A simple systematic method to proceed, without introducing the complication of inverted 
loops as in the standard spinor identity, can be found in [Q. First, use (43) of ||2^ (there 
is a factor of two missing on the left-hand side) to resolve any ip[a(3j5] into a sum of loop 
states depending on multiloops where a single loop contains at most three of the four loops 



a, P, 7, and 6. Then use (40) of |24| to rewrite any 'i/'[/3a7] in some preferred order ipla^-y] 
plus a sum of loop states depending only on multiloops where a single loop contains at most 
two of the three loops. Then our claim becomes 

ip[y U 7]tu U r]usi]st] - tph U 7?st U ijiuTjus] = 0, (74) 

which holds because 7 is inserted on the trivial loop, and the area derivative of a trivial loop 
dependence vanishes. Note that we used ( ^3|) for the area derivative at a kink several times. 

Given ([70[), the final observation is that in the commutator rj'' (t)'fj'^ {u) Abc{s) is anti- 
symmetric under exchange of t and u, while the remainder of the integral is symmetric under 
exchange of t and u. Therefore, the term in the integrand that does not reduce to a functional 
loop derivative vanishes under the integral. Hence we have shown that 

[Hss'iM),H,AN)] ^ J duoJaiv{u))gtiviu))j^y (75) 

The result represents one of the possible factor orderings of the operator version of D{ujag"'^), 
which appears in the classical Poisson bracket. 



5 Conclusion 

We have confirmed the result of that the constraint algebra of 3-1-1 quantum gravity 
in the loop representation formally closes without anomalies. There are two reasons for the 
comparative simplicity of our calculation. We were able to cast the Hamiltonian constraint 
into a form involving only functional loop derivatives instead of area derivatives, and we found 
a simple way to separate the rerouting operations from the other parts of the calculation. 

One direction for further research is to analyze the point-splitting regularization in more 
detail, in particular to analyze the next to leading order terms |Q. This appears to be 
necessary if one decides to take the point-splitting regularization seriously since one cannot 
remove the regulators in the definition of the constraints before computing their algebra with- 
out running into inconsistencies related to the background dependence of the regularization 
(also see the comment following (|6^). Turning this observation around, we avoid anoma- 
lous background dependent terms in the constraint algebra by postponing the removal of the 
regulators until the algebra is computed. 

A related point is that if we remove the regulator, the Hamiltonian constraint consists 
of discrete sums over kinks and intersections (no integrals involved) and integrals along the 
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loops for the acceleration terms ^, 20 1. The acceleration terms in a sense spoil the simple 
picture that the Hamiltonian constraint only acts on intersections, which are invariant under 
diffeomorphisms, although the background is present as an angle dependence, and the accel- 
eration terms depend on the background since they contain second derivatives of the loop, 
r/"(s). Sometimes one would like to argue away the acceleration terms, but notice that in 
order to obtain the integral on the right-hand side of the commutator of two Hamiltonians, 
@, and as is also apparent from section 4.3, it does not suffice to just consider the discrete 
sums corresponding to kinks and intersections. This may still be consistent with a diffeomor- 



phism invariant scheme like that of |11], in which acceleration terms do not appear, since, as 
we discussed in the introduction, (^) is no longer relevant. 

As a final remark, note that the rigorous framework based on diffeomorphism invariant 
measures |9|, is well adapted to the generators of diffeomorphisms, but has problems 

with the area derivatives appearing in the Hamiltonian constraint. Therefore it may be 



worthwhile to examine our form of the Hamiltonian constraint (30) in that setting. 

It is a pleasure to thank Abhay Ashtekar, John Baez, Roumen Borissov, Renate Loll and 
Carlo Rovelli for helpful discussions. 

References 

R. Jackiw, Quantum modifications to the Wheeler- De Witt equation, gr-qc/9506037| 



[2] 
[3] 
[4] 

[5] 

[6] 
[7] 
[8] 
[9] 



C. Rovelh and L. Smolin, Phys. Rev. Lett. 61 (1988) 1155; Nud. Phys. B331 (1990) 80 
A. Ashtekar, Phys. Rev. Lett. 57 (1986) 2244 

A. Ashtekar, Lectures on non-perturbative canonical gravity (World Scientific, Singapore, 
1991) 

B. Briigmann, Loop representations, in Canonical Gravity: From Classical to Quantum, 
eds. J. Ehlers and H. Friedrich, (Springer, 1994) 213 

M. Blencowe, Nucl. Phys. B341 (1990) 213 

R. Gambini, Phys. Lett. B255 (1991) 180 

B. Briigmann and J. Pullin, Nucl. Phys. B390 (1993) 399 

A. Ashtekar, J. Lewandowski, D. Marolf, J. Morao, and T. Thiemann, J. Math. Phys. 
36 (1995) 6456 



[10] Based on |||], the constraints can be obtained through the loop transform or directly 
from the loop operators. In the transform is used. In ||] the constraints are derived 
from the loop operators and shown to be equivalent to [^]. The very first derivation of 
the constraints in closed form in Q can be understood as a hybrid of the two methods, 
giving the same result (B. Briigmann, PhD Thesis Syracuse University, 1993). 

[11] C. Rovelh and L. Smohn, Phys. Rev. Lett. 72 (1994) 446 



18 



[12] A. Ashtekar and J. Lewandowski, Representation theory of analytic holonomy C* al- 
gebras, in Knots and Quantum Gravity, ed. J. Baez (Oxford University Press, Oxford, 
1994) 

[13] T. Thiemann, Reality conditions inducing transforms for quantum gauge field theory and 
quantum gravity, ^r-qc/9511057 ; A. Ashtekar, A generalized Wick transform for gravity, 
|gr-qc/951 10831 



[14] R. Gambini, A. Garat, and J. Pullin, The constraint algebra of quantum gravity in the 
loop representation, gr-qc/9404059 



[15] R. Gambini and J. Puhin, Loops, knots, gauge theories and quantum gravity (Cambridge 
University Press, in press) 

[16] J. Baez, Diffeomorphism invariant measures on the space of connections modulo gauge 
transformations, in Proc. of the Conf. on Quantum Topology, ed. D. Yetter (World 
Scientific, Singapore, 1994) 

[17] R. Loll, Nucl. Phys. B444 (1995) 619; Phys. Rev. Lett. 75 (1995) 3048 

[18] We thank A. Rendall for pointing out how a partition of unity allows one to extend the 
local integral of C^N to the manifold such that any test function M may be written as 
C^N for some v and A^. 

[19] B. Briigmann, R. Gambini, and J. Pullin, Phys. Rev. Lett. 68 (1992) 431. The tricky part 
is that one has to commute the limit which leads to diffeomorphism invariance with the 
limit defining the area derivative of the Hamiltonian constraint. A similar interchange 
of limits appears for example also in S. Major and L. Smolin, Quantum deformation of 
quantum gravity, |gr-qc/95 120201 

[20] R. Borissov, Regularization of the Hamiltonian constraint and the closure of the con- 
straint algebra, gr-qc/9411038| 

[21] N. Tsamis and R. Woodard, Phys. Rev. D36 (1987) 3691 

[22] J. Friedman and I. Jack, Phys. Rev. D37 (1987) 3495 

[23] B. Briigmann, R. Gambini, and J. Pullin, Nucl. Phys. B385 (1992) 587 

[24] B. Briigmann and J. Pullin, Nucl. Phys. B363 (1991) 221 



19 



